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(54) Recursive discrete fourier transformation 

(57) The invention is directed to an oxidizing com- 
position of agglutinated particles comprising an alkali 
metal pemnanganate, a binder/coating agent of an alkali 
metal silicate and/or alkali metal aluminate and a filler 
such as clay, diatomaceous earth, talc, ground silica or 



zeolites. The invention also is directed to an oxidizing 
kit which includes a pervious container and the oxidizing 
composition. The invention also is directed to a method 
of oxidizing waste materials with the kit in a substantially" 
aqueous moving fluid stream. 
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Description 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

[0001] The present invention relates to an arithmetic operation method for Fourier transformation or inverse Fourier 
transformation and more particularly to an arithmetic operation method for carrying out those transformations with a 
simple arithmetic operation processing for obtaining a result of the arithmetic operation by those transformations in a 
short time. 

2. Description of the Related Art 

[0002] Conventionally, the Fourier transformation method used in frequency analysis and the like of time-series data 
stream has been employed in high compression encoding of audio signal and moving picture signal as well as spectrum 
analysis of signals in audio signal processing field, moving picture data processing field for medical system and the 
like. Further, inverse Fourier transformation method and Fourier transformation method, which operate in complemen- 
tary way to each other, have been widely used as modulation technology and demodulation technology in communi- 
cation field. 

[0003] According to the Fourier transfonnation method, data stream sampled as digital amount is treated as N groups 
(N is an integer, for example, 1024) and a time interval, in which the N-data stream exist, is assumed as a window 
period while that window period is assumed to be basic frequency. Then, signal components of data stream existing 
to the window period are obtained as real number component and imaginary number component of a harmonic signal 
of the basic froquoncy. 

[0004] The inverse Fourier transfonnation method is a transfonnation method, which is complementary to the Fourier 
transformation method. 

[0005] Then, data stream supplied to be subjected to the Fourier transformation is discrete data stream sampled at 
every predetennined period and a method for conducting the Fourier transformation on the discrete data is called 
discrete Fourier transfonnation (DFT). The discrete Fourier transformation technology is employed as an analysis 
technology for control technology, in which, for example, a state of production process is obtained with discrete data, 
the quality of the process is kept optimum by analyzing the obtained data and the rate of acceptable products in man- 
ufactured products is improved. The inverse Fourier transfonriation method for discrete data is called inverse discrete 
Fourier transfonnation (IDFT) . which is employed for modulation of digital ground broadcasting, scheduled to be started 
from 2003 as well as generation of sound source in audio field, while its application field has been widening every year. 
[0006] The discrete Fourier transfomiation technology and the inverse discrete Fourier transformation technology 
which are employed in the above manner, are signal conversion technologies, which operate in complementary manner' 
In the discrete Fourier transfonnation of the fonner, supplied signals are sampled at every predetermined time interval 
and a voltage value obtained by the sampling is obtained as sampled data. When the data stream, which is a group 
of the obtained data, is N data obtained in a predetermined timet . x(t), x(t+1), x(t+2),..., x(t+N.2), x(t-hN-l), avaluex 
(k, t) of the discrete Fourier transfonnation obtained with respect to the N data is defined in a following equation 



50 



2 



cos 



2{n-t) 



^1 

N 



+ jc,(7i)sin 



2{n-t] 



jik' 

j7 



55 



j--^W- 



^ N 



+ X, 



COS 



2(„-0 



(1) 



for k = 0,1, .... N-1 



DOCID: <EP_:^1168192A2_L> 



2 



EP'l i68 192 A2 



[0007] As evident from this equation, the Fourier transformation is carried out on supplied data stream just such that 
its inherent basic function is convoluted in every point to be obtained and an operation for convoluting the basic function 
is conducted by a plurality of multiplication processings. 

[0008] When the multiplication processing is carried out with a special multiplying circuit or a digital signal processor 
5 (DSP), it has been known that a load upon hardware used for arithmetic operation In the multiplying circuit or the DSP 
is increased tremendously. 

[0009] As for the load on this multiplication processing, a transfonmation equation expressed with the aforementioned 
equation (1) needs multiplications of 4 X times, so that in case where N is 1024. for example, the multiplication 
frequency is about 4,200,000 times. Therefore, the load of the arithmetic operation is very large and further, if the 
10 number N of points to be handled as the data stream is increased, the multiplication processing frequency is increased 
by second power, which is unfavorable in terms of actual field. 

[0010] Thus, as the discrete Fourier transfonmation equation, fast Fourier transformation (FFT), whose operation 
efficiency is raised by deforming the matrix by paying attention to the regularity that the basic function is composed of 
a cyclic function, is employed actually. 
15 [0011] An arithmetic operation method called butterfly operation is employed for the FFT. The butterfly operation is 
so constructed that a simple integer, for example, 2 is defined as a basic number and supplied binary complex data Is 
subjected to complex operation comprised of a single addition, subtraction and multiplication each so as to output 
binary complex data. 

[0012] Therefore, the N-point FFT is comprised of log2 N-slep stage and (N/2)log2N butterfly operations and an 
20 operation result of the FFT can be obtained depending on the number of multiplications of (2N)log2N times. Conse- 
quently, the N-point FFT is used as a Fourier transformation method ensuring a high operating efficiency. 
[0013] The FFT or discrete Fourier transformation (DFT) conducted in the above manner is usually used, so that 
data stream, sampled successively at a predetermined time interval and then supplied, is subjected to Fourier trans- 
formation. In that Fourier transfomnation. the sampled data stream is stored in a memory circuit in succession tempe- 
rs rarily and when the number of data stored in the memory circuit temporarily reaches N, the Fourier transformation 
processing for that N data is started. 

[0014] During the Fourier transfomnation processing, data stored in the memory circuit temporarily is kept as it Is 
and when the arithmetic operation is terminated, new data is supplied to the memory circuit again so that it is stored 
therein. Then, after the predetemnined amount N of data is stored, the Fourier transfonnation processing is started 
30 again. 

[0015] However, because data is supplied without any interruption when the Fourier operation is executed also, a 
second memory circuit for accumulating N-data stream may be provided as another memory. Each time when N data 
are supplied, these memory circuits execute temporary storage processing and arithmetic operation processing alter- 
nately, so that the Fourier transformation is carried out on continuously supplied time-series data. 
35 [0016] In the Fourier transformation executed in such a manner, the arithmetic operation is carried out on N-data 
stream stored temporarily such that the data stream is treated as a block unit. Therefore, in such Fourier transformation 
processing, a delay of time corresponding to N-data stream is generated, so that a transformation result cannot be 
obtained at real time. 

[0017] Thus, although the data stream is supplied successively, the result of the Fourier analysis can be obtained 
40 only at N-samplings time interval. The delay time is not preferable for application of the Fourier analysis, thereby 
restricting its application field. 

[001 8] In order to secure a result of the Fourier transformation each time when sampled data is supplied, by executing 
the Fourier transformation processing on N-point data containing new data supplied successively to be sampled, to 
solve the above-described problem, N-point Fourier transformation processings need to be carried out within a sampling 

45 period. Even in case where FFT operation method developed for such fast arithmetic operation processing is employed, 
N-time operating speed is demanded to secure a Fourier operation result computed continuously from the FFT. In a 
normal case, il is difficult to conduct such super high speed FFT to supply an obtained operation result. 
[0019] Although in the Fourier transformation, generally, the FFT operation processing method is used, the FFT 
operation processing is carried out on N-point data stream, which is sampled with a clock signal whose sampling 

50 frequency is fs, quantitized and supplied, usually at a frequency interval of fs/N. 

[0020] In the arithmetic operation, data is fetched into other buffer memory even during the Fourier transformation 
processing, so that the Fourier analytic processing is carried out on time-serios data supplied during the operation 
also. In a period in which data stored in one buffer memory temporarily is being processed, data supplied to the other 
memory is fetched thereinto. When fetching of N data is temninated, those arithmetic operation and data fetching 

55 processing are switched so as to execute the Fourier operation processing. However, this method needs two sets of 
the buffer memories and PFT operation processing means, which is not a favorable method in economic viewpoint. 
[0021] Further, because this method carries out block processing in which supplied N time-series data are treated 
as a block, the result of the Fourier transformation on the N-point data fetched in is not outputted until N sampling clock 
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lime passes, and the analysis result obtained at that time is only a result of the Fourier transformation processing for 
every N samples. 

[0022] The result of the Fourier transfomnation on the newest N-point data containing data, which is sampled suc- 
cessively and supplied, cannot be outputted at real time. If it is intended to output the transformation result at real time, 
5 the Fourier transfomnation processing in each sampling period is needed. However, if the Fourier transformation is 
carried out continuously at every sampling time interval, the amount of arithmetic operation per unit time becomes 
tremendous, which is not a realistic method. 

[0023] On the other hand, as a method for obtaining a result of continuous Fouriertransformation processing, "Fourier 
transformation device and Fouriertransformation method" has been proposed in Japanese Patent Aoplication Laid- 
10 Open No, H1 -59454. 

[0024] This publication describes a method for Fourier-transforming a vibration waveform value sampled and sup- 
plied. According to the same method, a difference between a newly supplied vibration waveform value and an old 
vibration wavefomi value already supplied and used for Fourier transfomiation processing is obtained and a newer 
complex vibration value than the old complex vibration value provided by the arithmetic operation is obtained each 
f5 time when a sampling value of the vibration wavefomi is supplied. 

[0025] Although frequency analysis may be considered as an application using the continuous Fourier transfonna- 
tion, if such application is considered, often any frequency band is analyzed at any resolution and then, that application 
is executed depending on the result. 

[0026] However, the aforemenlioned publication does not consider those points. 
20 [0027] Input data stream described in the same publication refers to only real number. Although if the data stream 
to be sampled is data stream which can be expressed in time series, it is permitted that the data stream is limited to 
the real number, imaginary number input signal as well as real number input signal may be sometimes used in Fourier 
operation for use in communication field, and therefore, supplied complex signals have to be subjected to transformation 
processing. 

[0028] However, the conventional fast Fourier transformation method, in which continuously supplied data is sub- 
jected to Fourier transformation to secure an operation result, is not capable of carrying out complex Fourier transfor- 
mation on data stream supplied in the fomn of complex number. 
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SUMMARY OF THE INVENTION 



[O029] Accordingly, an object of the present invention is to provide a continuous Fourier transformation method suit- 
able for transfomning any frequency band desired to be subjected to Fourier analysis at a desired resolution. 
[0030] Another object of the present invention is to provide a continuous Fourier transfomnation method enabling 
Founer analysis on data stream supplied in the form of complex number and suitable fortransfomiing any frequency 
35 band desired to be subjected to Fourier analysis at a desired resolution. 

[0031] To achieve the above object, there is provided a recursive discrete Fourier transfomiation method wherein 

data values x(t), x(t+1 ), x(t+2), x{t+3) x(t+N-1 ), x(t+N) sampled at times t, t+1 , t+2, t+3. .... t+N-1 . t+N (N is a positive 

integer which is 1 or more) each having an equal interval are supplied and with such N data values supplied since time 
t as a data stream, a frequency component, which is degree k (k is 0 or a positive integer smaller than N) obtained by 
carrying out complex Fourier transformation on the data stream, is obtained such that a real part Xr(k. t) thereof and 
an imaginary part Xi(k, t) thereof are complex Fourier coefficients, the method comprising: a first step of storing the 
data stream x(t). x(t+1 ). x(t+2). x(t+3), x(t+N-1 ) supplied since time t at time t+N-1 into a first memory means tem- 
poranly; a second step of obtaining the complex Fourier coefficients X,(k, t) and Xi(k, t) of the data stream stored 
temporanly at the first step; a third step of storing the complex Fourier coefficients X,{K t) and Xi(k, t) obtained at the 
second step Into the second memory means temporarily; and a fourth step of by using data value x(t+N) supplied at 
time t+N, data value x(t) stored in the first memory means temporarily and the complex Fourier coefficients X,(k, t), Xj 
(k.t) stored in the second memory means temporarily, obtaining complex Fourier coefficients X,(k, t+1 ) and Xi(k, t+1 ) 
to the data stream supplied since time t+1 with respect to a positive constant value A for giving an amplitude value 
corresponding to a difference between the x(t+N) and the x(t) according to a following equation: 



xAk,t +1)= f Ar)-4)]|xcos|2^j j:,(;t,r)sm |2^j 
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[0032] According to the present invention, data stream comprised of supplied N sannples is subjected to connplex 
Fourier transfonmation and then, the oldest sample value used for complex Fourier transfonnation is subtracted from 
the value of a new sample to be supplied next so as to obtain a subtraction value. By deleting the oldest sample at the 

10 same time, discrete Fourier operation result on new N samples can be obtained based on that subtraction value and 
the result of the preceding complex Fourier operation. Thus, different from the conventional Fourier operation that Is 
carried out after N sample data is supplied, a Fourier operation result can be obtained quickly in a single sample period. 
[0033] Further, to achieve the above object, there is provided a recursive discrete Fourier transfonnation method 
wherein with sampling frequency as fs, data values x(t), x(t+1), x(t+2), x(t+3), .... x(t+N-1), x(t+N) sampled at times t, 

15 t+1 , t4-2, t+3, .... t+N-l . t+N (N is a positive integer which is 1 or more) each having an equal interval are supplied and 
with such N data values supplied since time t as a data stream, a frequency interval given with a iriinirnum frequency 
' f 1 and a maximum frequency f2 to the data stream is regarded as a measuring frequency interval, airequency Interval 
'Obtained "by dividing that measuring frequency interval by the N is assumed to be an analysis frequency interval and 
" a result of the frequency analysis obtained by carrying out complex Fourier transformation al every analysis frequency 

20 interval is obtained as a frequency component which is k (k is 0 or a positive integer smaller than N) times the frequency 
interval, such that a real part Xr(k, t) thereof and an imaginary part Xj{k, t) thereof are complex Fourier coefficients, the 
method comprising: a first step of storing the data stream x(t), x(t+1), x(t+2), x(t+3), x(t+N-1) supplied since time t 
at time t+N-l into a first memory means temporarily; a second step of obtaining the complex Fourier coefficients X^{k, 
t) and Xj(k, t) of the data stream stored temporarily at the first step; a third step of storing the complex Fourier coefficients 

25 Xr(k, t) and Xj(k, t) obtained at the second step into the second memory means temporarily; and a fourth step of by 
using data value x(t+N) supplied at time t+N, data value x(t) stored in the first memory means temporarily and the 
complex Fourier coefficients X^(k, t). Xj(k,t) stored in the second memory means temporarily, obtaining complex Fourier 
coefficients X^(k, t+1) and Xj(k, t+1) within the frequency interval with the minimum frequency f1 and the maximum 
frequency f2 to the data stream supplied since time t+1 with respect to a positive constant A for giving an amplitude 

30 corresponding to a difference between the x(t4-N) and the x(t) according to following equations: 
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[0034] According to the present invention, data stream comprised of supplied N samples is subjected to complex 
Fourier transformation and then, the oldest sample value used for complex Fourier transfonnation is subtracted from 
the value of a new sample to be supplied next so as to obtain a subtraction value. By deleting the oldest sample at the 
same time, discrete Fourier operation result on new N samples can be obtained based on that subtraction value and 
the result of the preceding complex Fourier operation. At this time, a frequency band for analysis for the discrete Fourier 
operation is set up with f1 and f2 preliminarily and the Fourier operation is carried out with a basic frequency whose 
resolution is set to N. Thus , different from a conventional case in which the frequency band and resolution are auto- 
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matically set up when the Fourier operation is carried out after the N sample data is supplied, a result of the Fourier 
operation within any frequency band at any resolution can be obtained quickly in a single sampling period. 
[0035] According to a preferred embodiment of the present invention, inverse discrete Fourier transfomnation data 
subjected to inverse discrete Fourier transformation using a positive constant a is supplied so as to obtain a complex 
Fourier transformation coefficient to the supplied inverse discrete Fourier transformation data, the discrete Fourier 
transformation being carried out using the positive constant A corresponding to the positive constant a. 
[0036] According to this embodiment, in addition to the above described effect, the operation can be carried out by 
FFT in which the positive constant A is set up by selecting appropriately 1 , N or a square root of N. This is employed 
in. for example, communication apparatus. In case where information signal to be transmitted is subjected to inverse 
fast Fourier transformation (IFFT) processing to be converted to time-series signal and that time-series signal is trans- 
mitted, in Fourier operation method for obtaining information signal by supplying that transmitted signal to the FFT for 
executing the FFT operation, the constant value A corresponding to the constant value a, which is, for example, 1 , N 
or the square root of N used for the IFFT operation, is employed in the FFT operation processing circuit, so as to 
construct the FFT which operates in a complementary way with the IFFT to decode the information signal. By using 
the constant value A corresponding to the value a used for the IFFT operation for a signal supplied from other mating 
system, system having an excellent characteristic can be established. 

[0037] According to a preferred embodiment of the present invention, k indicating a degree of the complex Fourier 
coefficient is set to a desired value and Fourier transformation to the set value k is repeated at every time having an 
equal interval. 

[0038] According to this embodiment, in addition to the above described effect, particularly by setting up the value 
of k to a specific value, a result of the Fourier analysis on a specific frequency determined with the value k can be 
computed quickly. Thus, if there is a system controlled according to the result of the Fourier transformation on that 
specific frequency, the control of that system can be earned out at a minimunn delay time, so that response characteristic 
of the entire system can be improved. Thus, system having a better quality can be constructed. 
[0039] Further, to achieve the above object, there is provided a recursive discrete Fourier transformation method 

wherein data values x(t), x(t+1), x(t+2), x(t+3),..., x{t+N-1), x(t+N) sampled at times t, t+1 , t+2, t+3 t+N-1 , t+N (N 

is a positive integer which is 1 or more) each having an equal interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which is degree k (k Is 0 or a positive integer smaller than N) 
obtained by carrying out complex Fourier transformation on the data stream, is obtained such that a real part X^(k, t) 
thereof and an imaginary part Xj(k, t) thereof are complex Fourier coefficients, the method comprising: a first step of 
storing the data stream x(t), x(t+1 ), x(t+2), x(t+3) , x(t+N-1 ) supplied since time t at time t+N-1 into a first memory 
means temporarily; a second step of obtaining the complex Fourier coefficients X^{k, t) and Xi(k, t) of the data stream 
stored temporarily at the first step; a third step of storing the complex Fourier coefficients X^{k, t) and Xj(k, t) obtained 
at the second step into the second memory means temporarily; and a fourth step of by using the complex Fourier 
coefficients X^{k, t) - jXj(k, t) stored in the second memory means temporarily, obtaining complex Fourier coefficients 
Xr(k, t+1 ) - jXj(k, t+1 ) to the data stream supplied since time t+1 based on a transfer function expressed in a following 
equation. 
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where A is a positive constant for providing [x(t+N) - x(t)] with an amplitude. 

[0040] According to the present invention, data stream comprised of supplied N samples is subjected to complex 
Fourier transformation and then, the oldest sample value used for complex Fourier transformation is subtracted from 
the value of a new sample to be supplied next so as to obtain a subtraction value. By deleting the oldest sample at the 
same time, discrete Fourier coefficients to new N samples can be obtained by operation with a predetemnined transfer 
function based on that subtraction value and the preceding complex Fourier operation result. Thus, different from the 
conventional Fourier operation that is carried out after N sample data is supplied, a Fourier operation result can be 
obtained quickly in a single sample period. 

[0041] Further, to achieve the above object, there is provided a recursive discrete Fourier transformation method 
wherein with sampling frequency as fs, data values x(t). x(t+1), x(t+2), x(t+3), x(t+N-1), x(t+N) sampled at times t, 
t+1 , t+2, t+3, t+N-1 , t+N (N is a positive integer which is 1 or more) each having an equal interval are supplied and 
with such N data values supplied since time t as a data stream, a frequency interval given with a minimum frequency 
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f 1 and a maximum frequency f2 to that data stream is regarded as a measuring frequency interval, a frequency interval 
obtained by dividing that measuring frequency interval by the N is assumed to be an analysis frequency interval and 
a result of the frequency analysis provided by carrying out complex Fourier transformation at every analysis frequency 
interval is obtained as a frequency component which is k (k is 0 or a positive integer smaller than N) times the frequency 

5 interval, such that a real part Xr(k, t) thereof and an imaginary part Xj(k, t) thereof are complex Fourier coefficients, the 
method comprising: a first step of storing the data stream x(t), x(t+1), x(t+2), x(t+3), .... x(t+N-1) supplied since time t 
at time t+N-1 into a first memory means temporarily; a second step of obtaining the complex Fourier coefficients X^(k, 
t) and Xi(k, t) of the data stream stored temporarily at the first step; a third step of storing the complex Fourier coefficients 
Xr(k, t) and Xj(k, t) obtained at the second step into the second memory means temporarily; and a fourth step of by 

10 using the complex Fourier coefficient X^(K t) - jXj(k, t) stored in the second memory means temporarily, obtaining 
complex Fourier coefficients X^(K t+1) - jXj (k, t+1) at the frequency interval given by the minimum frequency f1 and 
the maximum frequency f2 to the data stream supplied since time t+1 based on a transfer function expressed in a 
following equation: 

ym^^M^:!im0: ^^ „.:..■ 

20 

where A is a positive constant for providing [x(t+N) - x(t)) with an amplitude value, and 
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[0042] According to the present invention, data stream comprised of supplied N samples is subjected to complex 

30 Fourier transformation and then, the oldest sample value used for complex Fourier transfomnation is subtracted from 
the value of a new sample to be supplied next so as to obtain a subtraction value. By deleting the oldest sample at the 
same time, discrete Fourier operation result on new N samples can be obtained based on that subtraction value and 
the result of preceding the complex Fourier operation . At this time, a frequency band for analysis for the discrete Fourier 
operation is set up with f1 and f2 preliminarily and that Fourier operation is carried out with a basic frequency whose 

35 resolution is set to N. Thus, different from a conventional case in which the frequency band and resolution are auto- 
matically set up when the Fourier operation is carried out after the N sample data is supplied, a result of the Fourier 
operation within any frequency band at any resolution can be obtained quickly in a single sampling period. 
[0043] According to a preferred embodiment of the present invention , the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) and the x(t) is capable of being set selectively with an 

40 inverse number of square root of 1 , N or 1/N. 

[0044] According to this embodiment, in addition to the above described effect, the operation can be carried out by 
FFT in which the positive constant A is set up by selecting appropriately 1 , 1/N or an inverse number of a square root 
of N. This is employed in, for example, communication apparatus. In case where infonnation signal to be transmitted 
is subjected to inverse fast Fourier transfonnation (I FFT) processing to be converted to time-series signal and that 

45 time-series signal is transmitted, in Fourier operation method for obtaining information signal by supplying that trans- 
mitted signal to the FFT for executing the FFT operation, the constant value A con-esponding to the constant value a, 
which is, for example. 1. 1/N or the inverse number of the square root of N used for Ihe IFFT operation, is employed 
in the FFT operation processing circuit, so as to construct the FFT which operates in a complementary way with the 
IFFT to decode the information signal. By using the constant value A corresponding to the value a used for the IFFT 

50 operation for a signal supplied from other mating system, system having an excellent characteristic can be established. 
[0045] According to a preferred embodiment of the present invention, k indicating the degree of the complex Fourier 
coefficient is set to a desired value and the Fourier transformation to the set value k is repeated at every time having 
an equal interval. 

[0046] According to this embodiment, in addition to the above described effect, particularly by setting up the value 
55 of k to a specific value, a result of the Fourier analysis on a specific frequency determined with the value k can be 
computed in a short time. Thus, if there is a system controlled according to the result of the Fourier transfomnation on 
that specific frequency, the control of that system can be carried out at a minimum delay time, so that response char- 
acteristic of the entire system can be improved. Thus, system having a better quality can be constructed. 
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[0047] Further, to achieve the above object, there is provided a recursive discrete Fourier transformation method for 
obtaining complex Fourier coefficients, wherein data values x,(t)+jXi(t), x,(t+1 ) + jXi(t+1 ), x,(X^2) + jXj(t+2), Xr(t+3) + jXj 
{t+3). .... x,(t+N-1) +iXi(t+N-1), x,(t+N) jXi(t+N) sampled at times t, t+1 , t+2,t-t-3, t+N-1 , t+N (N is a positive integer 
which is 1 or more) each having an equal interval are supplied and with such N data values supplied since time t as a 
data stream, a frequency component, which is degree k (k is 0 or a positive integer smaller than N) obtained by carrying 
out complex Fourier transformation on the data stream, is obtained such that a real part X,{K t) thereof and an imaginary 
part Xj(k t) thereof are complex Fourier coefficients, the method comprising: a first step of storing the data stream x^ 
(t)+jXi(t), x,(t+1) +jXj(t+1). x,(t+2) + jXj(t+2). x,(t4-3) + jXi(t+3). x,{t+N-1) + jXi(t+N-1) supplied since time t at time t+N- 
1 into a first memory means temporarily; a second step of obtaining the complex Fourier coefficients X,{K t) and Xj(k, 
t) of the data stream stored temporarily at the first step; a third step of storing the complex Fourier coefficients X,(k, t) 
and Xi(k, t) obtained at the second step into the second memory means temporarily; and a fourth step of by using the 
data value x^a+N) + jXj(t+N) supplied at time t+N, the data value x,(t) +JXj(t) stored in the first memory means temporarily 
and the complex Fourier coefficient X,(k, t) and Xi(k, t) stored in the second memory means temporarily and used 
recursively, obtaining complex Fourier coefficients t+1 ) and Xi(k, t+1 ) to the data stream supplied since time t+1 
with respect to a positive constant A for giving an amplitude value corresponding to the difference between the Xr(t+N) 
and the x^{X) according to a following equation: * 
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[0048] According to the present invention, data stream comprised of supplied N complex data is subjected to Fourier 
transformation, so that a Fourier operation result to new N complex data can be obtained based on the preceding 
Fourier transfomiation result, data value of new complex data to be supplied next and the oldest complex data used 
for executing the Fourier transformation. Therefore, different from the conventional case where the Fourier operation 
is carried out after the N sample data is supplied, the Fourier operation result can be obtained quickly in the period of 
a single complex data. 

[0049] Further, to achieve the above object, there is provided a recursive discrete Fourier transfomnation method for 
obtaining complex Fourier coefficient wherein with sampling frequency as fs, data values Xr(t)+jXj(t), Xr(t+1) + jXj(t+1). 

Xr(t+2) + jXi(t+2). Xr(t+3) + JXj(t+3) x,(t+N-1) + jXi(t+N-1). x,(t+N) + jXj(t+N) sampled at times t, t+1 , t+2, t+3. .... t+N- 

l' l+N (N is a positive integer which is 1 or more) each having an equal interval are supplied and with such N complex 
data values supplied since time t as a data stream, a frequency interval given with a minimum frequency f1 and a 
maximum frequency f2 to that data stream is regarded as a specified frequency interval, a frequency interval obtained 
by dividing the specified frequency internal by the N is assumed to be a minimum frequency interval and a result of 
the Fourier transformation provided by carrying out complex Fouriertransfomnation at every minimum frequency interval 
is obtained as a frequency component which is k (k is 0 or a positive integer smaller than N) times the minimum 
frequency interval, such that a real part X,{k, t) thereof and an imaginary part Xj(k, t) thereof are complex Fourier 
coefficients, the method comprising: a first step of storing the data stream x,(\) + jXi{t), x^O+l) + jx,(t+1), x,{t+2) + jx, 
(t+2), Xr(t+3) +jXj(t+3), .... x,(t+N-1) + jXi{t+N-1). Xr(t+N) + jXi(t+N) supplied since time tat time t+N-1 into a first memory 
means temporarily: a second step of obtaining the complex Fourier coefficients X,(K t) and Xj(k, t) of the data stream 
stored temporarily at the first step; a third step of storing the complex Fourier coefficients X,(k, t) and Xi(k, t) obtained 
at the second step into the second memory means temporarily; and a fourth step of by using the data value x,(t+N) + 
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jXj(t+N) supplied at time t+N, the data value x^(\) + jXf(t) stored temporarily in the first memory means and complex 
Fourier coefficients X^iK t) and Xj(k, t) stored temporarily in the second memory means and used recursively, obtaining 
complex Fourier coefficient X^(k, t+1) and Xi(k. t+1) in the specified frequency interval given by the minimum frequency 
f 1 and the maximum frequency f2 to the data stream supplied since time t+1 with respect to a positive constant A for 
5 giving an amplitude value corresponding to a difference between the x(t+N) and the x(t) according to. a following equa- 
tion: 




30 [0050] According to the present invention, data stream comprised of supplied N complex data is subjected to Fourier 
transformation and then, when it is intended to obtain a result of the Fourier operation to new N complex data on the 
basis of the result of the preceding Fourier transfonnation, data value of new complex data to be supplied next and 
the oldest complex data used for Fourier, the frequency band to be Fourier-transformed is set up with f 1 and f2 pre- 
liminarily and further, the Fourier operation is can-ied out using the basic frequency set with the resolution N. Therefore, 

35 different from the conventional case where the frequency band and the resolution are automatically determined when 
the N sample data are supplied for the Fourier operation, a Fourier operation result to any frequency band and resolution 
can be obtained quickly and easily in the period of a single sample. 

[0051] According to a preferred embodiment of the present invention, by using data value y^{X) + jyj(t) supplied at 
time t, data value yr(t+N) + jyj(t+N) supplied at time t+N and complex inverse discrete Fourier coefficients Yr(k, t) and 
40 Yi(k, t) with the real part and the imaginary part obtained with respect to N data values )/^(X) + jyj(t), yr(t+1) + jyj(t+1). .... 
yr(t+N-1 ) + jyj(t+N-1 ) supplied from time t to time t+N-1 , a real part Y^Ck. t+1 ) and an imaginary part Yj(k, t+1 ) of each 
complex inverse discrete Fourier coefficient of N data values supplied since time t+1 are obtained with respect to a 
positive constant value B for giving an amplitude corresponding to a difference between the yr(t+N) and the y^(\) as 
inverse discrete Fourier transfonnation data according to following equations, 

45 
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and after the obtained inverse discrete Fourier transfomnation data is supplied, discrete Fourier transfomnation on the 
supplied inverse discrete Fourier transfomnation data is carried out, the recursive discrete Fourier transformation con- 
ducting discrete Fourier transfomnation by using the positive constant value A corresponding to the positive constant 
value B. 

[0052] According to this embodinnent, in addition to the above described effect, the Fourier transfomnation can be 
carried out such that the positive constant A is set up to 1 . N or a square root of N. This is employed in, for example, 
communication apparatus. In case where information signal to be transmitted is subjected to inverse fast Fourier trans- 
formation (IFFT) processing to be converted to time-series signal and that time-series signal is transmitted, in the 
Fourier operation method for obtaining information signal by supplying that transmitted signal to the FFT for executing 
the FFT operation, the constant value A corresponding to the constant value B, which is, for example, 1, N or the 
square root of N used for the IFFT operation, is employed in the FFT operation processing circuit, so as to construct 
the FFT which operates in a complementary way with the IFFT. By using the constant value A corresponding to the 
value B used for the IFFT operation for a signal supplied from other mating system, system containing Fourier trans- 
formation processing having an oxcelient characteristic can be established. 

[0053] According to a preferred embodiment of the present invention, k indicating the degree of the complex Fourier 
coefficient is set to a desired value and the Fourier transfomnation to the set value k is repeated at every time having 
an equal interval. 

[0054] According to this embodiment, in addition to the above described effect, particulariy by setting up the value 
of k to a specific value, the Fourier analysis on a specific frequency determined with the value k can be carried out in 
a short time. Thus, if there is a system controlled according to the result of the Fourier transfomnation on that specific 
frequency, the control of that system can be carried out at a minimum delay time, so that response characteristic of 
the entire system can be improved. Thus, system whose response characteristic is a short time, having a better quality, 
can be constructed. 

[0055] Further, to achieve the above object, there is provided a recursive inverse discrete Fourier transfomnation 
method for obtaining complex inverse Fourier coefficient wherein complex data values x^(X) + jXj(t), x,(t+1 ) + jXi(t+1 ), x^ 

(t+2) + jXj(t+2), Xr(t+3) + jXj{t-f3) Xr(t+-N-1 ) + jXj(t+N-1 ), x,{t+H) + jXi(t+N) sampled at times t, t+1 . t+2, t+3. t+N-1 , 

t+N (N is a positive integer which is 1 or more) each having an equal interval are supplied and with such N data values 
supplied since time t as a data stream, an inverse complex Fourier transfomnation component, which is degree k (k is 
0 or a positive integer smaller than N) obtained by carrying out inverse complex Fourier transfomnation on the data 
stream, is obtained such that a real part t) thereof and an imaginary part Xj(k. t) thereof are inverse complex 
Fourier coefficients, the method comprising: a first step of storing the data stream Xr(t)+jXj(t), Xr(t+1) + jXj(t+1), x^(X+2) 
+ jXj(t+2), \{t+3) + jXj(t+3), Xr(t-i-N-l) + jXj(t+N-1) supplied since time t at time t+N-1 into a first memory means 
temporarily; a second step of obtaining the inverse complex Fourier coefficients Xr(k, t) and Xj(k, t) of the data stream 
stored temporarily at the first step; a third step of storing the inverse complex Fourier coefficients X^(K t) and Xi(k, t) 
obtained at the second step into the second memory means temporarily; and a fourth step of by using the data value 
X|.(t+N) -I- jXj(t+N) supplied at time t+N, the data value x,(t) + jXi(t) stored in the first memory means temporarily and the 
inverse complex Fourier coefficient X^(k, t) and Xi{k, t) stored in the second memory means temporarily obtaining 
inverse complex Fourier coefficients X^(k, t+1) and Xj(k. t+1) to the data stream supplied since time t+1 with respect 
to a positive constant B for giving an amplitude value corresponding to the difference between the Xr(t+N) and the x^ 
(t) according to a following equation: 
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[0056] According to the present invention, data stream comprised of supplied N complex data is subjected to inverse 
Fourier transformation, so that an inverse Fourier operation result to new N complex data can be obtained based on 
the preceding inverse Fourier transformation result, data value of new complex data to be supplied next and the oldest 
complex data used for executing the inverse Fourier transformation. Thcreforo, different from the conventional case 
25 where the inverse Fourier operation Is carried out after the N sample data is supplied, the inverse Fourier operation 
result can be obtained quickly in the period of a single complex data. 

[0057] Further, to achieve the above object, there is provided a recursive inverse discrete Fourier transformation 
method for obtaining complex inverse Fourier coefficients wherein with sampling frequency as fs, data values x^(t) + 
jXj(t). x,(t+1) + jXi(t+1), x^(t+2) + jXj(t+2), x^(t+3) + jXj(t+3), .... x,(t+N-1) + jXj(t+N-1), x^Ct+N) + jXj(t+N) sampled at times 

30 t, t+1, t+2, t+3, .... t+N-1 , t+N (N is a positive integer which is 1 or more) each having an equal interval are supplied 
and with such N data values supplied since time t as a data stream, a frequency interval given with a minimum frequency 
f 1 and a maximum frequency f2 to the data stream is regarded as a specified frequency interval, a frequency interval 
obtained by dividing that specified frequency interval by the N is assumed to be a minimum frequency interval and a 
result of the inverse Fourier transformation provided by carrying out inverse complex Fourier transfomnatlon at every 

35 minimum frequency interval is obtained as a frequency component which is k (k is 0 or a positive integer smaller than 
N) times the minimum frequency interval, in the form of a real part Xr(k, t) and an imaginary part Xj(k, t), the method 

comprising: a first step of storing the data stream Xr(t)+-jXi(t) , Xr(t+1) + jXj(t+1), Xr(t+2) + jXi(t+2): Xr(t+3) + jXj(t+3) x^ 

(t+N-1) + jXj(t+N-1), Xr(t+N) + jXj(t+N) supplied since time t at time t+N-1 into a first memory means temporarily; a 
second step of obtaining the inverse Fourier coefficients Xj.(k, t) and Xj(k, t) of the data stream stored temporarily at 

40 the first step; a third step of storing the inverse Fourier coefficients Xr(k, t) and Xi(k, t) obtained at the second step into 
the second memory means temporarily; and a fourth step of by using the data value Xr(t+N) + jXj(t+N) supplied at time 
t+N, the data value x^{Vj + jXj(t) stored temporarily in the first memory means and preceding inverse Fourier coefficients 
Xr(k, t) and Xj(k, t) stored temporarily in the second memory means, obtaining inverse Fourier coefficient X^lk, t+1 ) and 
Xj(k, t+1) in the specified frequency interval given by the minimum frequency f1 and the maximum frequency t2 to the 

45 data stream supplied since time t+1 with respect to a positive constant B for giving an amplitude value corresponding 
to a difference between the x(t+N) and the x(t) according to following equations: 
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[0058] According to the present invention, data stream comprised of supplied N complex data is subjected to inverse 
Fourier transformation and then, when it is intended to obtain a result of the inverse Fourier operation to new N complex 
data on the basis of the result of the preceding inverse Fourier transfomriation, data value of new complex data to be 
supplied next and the oldest complex data used for inverse Fourier, the frequency band to be inverse-Fourier-trans- 
15 formed is set up with f1 and f2 preliminarily and further, the inverse Fourier operation is carried out using the basic 
frequency set with the resolution N. Therefore, different from the conventional case where the frequency band and the 
resolution are automatically determined when the N sample data are supplied for the inverse Fourier operation, an 
inverse Fourier operation result to- any frequency band and resolution can be obtained quickly and easily in the period 
of a single sample. 

20 [0059] According to a preferred embodiment of the present invention, k indicating the degree of the complex Fourier 
coefficient is set to a desired value and the inverse Fourier transformation to the set value k is repeated at every time 
having an equal Interval. 

[0060] According to this embodiment, in addition to the above described effect, particularly by setting up the value 
of k to a specific value, the inverse Fourier analysis on a specific frequency determined with the value k can be carried 
25 out in a short time. Thus, if there is a system for generating a signal according to the result of the inverse Fourier 
transformation on that specific frequency, an operation of that system can be executed at a minimum delay time, so 
that response characteristic of the entire system can be improved. Thus, because the response characteristic of the 
entire system can be improved, system whose response characteristic is shorter, having an excellent quality, can be 
constructed. 

30 [0061] The nature, principle and utility of the invention will become more apparent from the following detailed de- 
scription when read in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

35 [0062] In the accompanying drawings: 

Fig.1 is a diagram showing schematically a recursive discrete Fourier transformation device according to an em- 
bodiment of the present invention; - 

Fig.2 Is a diagram showing schematically a relation between data value obtained by sampling the waveform of 
40 signal supplied to the recursive discrete Fourier transformation device of the embodiment of the present invention 

and a corresponding DFT operation; 

Fig.3 is a diagram showing the structure of the recursive discrete Fourier transfomriation, which carries out arith- 
metic operation using a transfer function according to the embodiment of the present invention; 
F(g.4 is a diagram showing an analysis frequency band in the recursive discrete Fourier transformation device of 
45 the embodiment of the present invention and a relation of the frequency relating to setting of resolution; 

Fig. 5 is a diagram showing the structure for an analysis frequency band and the setting of the resolution in the 
recursive discrete Fourier transfomriation device which carries out arithmetic operation using a transfer function of 
the embodiment of the present invention; and 

Fig. 6 is a diagram showing a schematic structure of the recursive inverse discrete Fourier transformation device 
50 of the embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0063] Hereinafter, the preferred embodiment of the recursive discrete Fourier transformation method and the recur- 
55 sive inverse discrete Fourier transformation method of the present invention will be described with reference to the 
accompanying drawings. 

[0064] Fig.1 shows an embodiment of a recursive discrete Fourier transformation device applying the recursive dis- 
crete Fourier transformation method, which will be described below with reference to the same Figure. 
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[0065] The recursive discrete Fourier trRnsformation device comprises a data updating portion 1 to which data sam- 
pled in a predetermined time interval Is supplied and which stores temporarily supplied new N (N: positive integer) data 
therein, a basic frequency setting portion 2 in which the basic frequency for carrying out the discrete Fourier transfor- 
mation Is set up, a recursive DFT operating portion 3 for carrying out recursive DFT operation, and a memory 4 for 
5 storing temporarily computed data. 

[0066] Next, an operation of the recursive discrete Fourier transformation device having such a structure will be 
described. 

[0067] First, supplied data is sampled with a sampling circuit (not shown) every predetermined time interval and 
then, the sampled and quantltized discrete data is supplied to the data updating portion 1 . 

10 [0068] The sampling circuit samples data values supplied at every time of a predetennined interval, 1. 1+1 , X+2, t+3, .. 
t+N-1, t+N (N: natural number) and generates data streanr^ x(t), x(t+1). x(t+2). x{t+3), .... x(t+N-1), x(t-i-N) by converting 
the data values supplied at the respective times to sampling values corresponding to the respective times. 
[0069] The operation of the sampling circuit is the same as that of an A/D converter which converts supplied analog 
signal to digital signal, converting an analog signal voltage value supplied at a predetermined time interval, given in 

15 the form of an inverse number of a sampling frequency, to a digital signal value. The converted digital signal value is 
a voltage having an analogous relation with an analog voltage, which gives a pulse amplitude modulation signal or that 
voltage value is expressed with a binary digital value. 

[0070] In this manner, sampling datax(t) sampled at time t, x(t+1) sampled at time t+1, x(t+N-1) sampled att+N- 
1 , x(l+N) sampled at lime l+N, .... are supplied Lo the data updating portion 1 . 
20 [0071] The data-updating portion 1 updates the newest N data (N: positive integer) of the supplied data, and stores 
them temporarily. 

[0072] That is, data supply is started from x(t) and when data x(t+1 ), x(t+2), x(t-i-3) are supplied, all the supplied data 
x(t), x(t+1), x(t+2), x(t+3) are stored temporarily. This temporary storage operation is continued until data x(t-t-N-l) is 
entered and when data x(tH N-1 ) is entered, the total data number becomes N, so that data region of the data updating 

25 ' p6rtion 1 becomes full. 

[0073] When next data x(t+N) is supplied with such a condition, the total data number becomes N+1 . Thus, the data 
updating portion 1 subtracts x(t) from x{t+N) and then supplies data obtained by that subtraction to the recursive DFT 
operating portion 3 while removing the oldest data x(t) from the memory portion. Consequently, the data updating 
portion 1 temporarily stores N data x(t+1), x(t+2), x(t+3), .... x(t+N-1), x(t-i-N). 

30 [0074] Likewise, when next data x(t+N+1 ) is supplied, the data updating portion 1 obtains x(t-t-N+1)-x(t+1 ) and sup- 
plies it to the recursive DFT operating portion 3 while removing x(t+1) from the memory. Consequently, the newest N 
data of supplied data is always stored in the data-updating portion 1 temporarily. 

[0075] The N data stored temporarily is supplied to the recursive DFT operating portion 3 and a last FFT operation 
result stored in the memory 4 temporarily is supplied to the recursive DFT operating portion 3 as recursive data. Then, 
35 the recursive DFT operating portion 3 carries out recursive discrete Fourier transfomnation according to frequency 
resolution information set up by the basic frequency setting portion 2 in a manner described later and outputs a result 
of Its arithmetic operation. 

[0076] Next, this recursive discrete Fourier transfomnation operation method will de described in detail with reference 
to a conventional Fourier operation as well. 
40 [0077] Fig.2 explains schematically a relation between data value obtained by sampling the wavefomn of supplied 
signal In a sampling period ts and corresponding DFT operation. 

[0078] The same Figure shows a set of N real number data values x(t), x(t+1), x(t+2), x(t+3), .... x(t+N-1) sampled 
in a predetermined sampling period from time t and a set of N real number data values x{t+1), x(t+2), x(t+3), x(t+N- 
1 ), x(t+N) sampled from time t+1 . 
45 [0079] Then, the value X(k, t) of discrete Fourier transfomnation obtained with respect to data stream which is N real 
number data values sampled since time t is defined in the following expression. 
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[0080] Transformation of real part X^(K t) and imaginary part Xj(k, t) of complex Fourier transformation of sampling 
data stream x(t) ~ x(t+N-1 ) supplied since arbitrary time t is defined as follows. 
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[0081] The transformation equation for the sampling data supplied since timet is defined In the above manner. Next, 
a transformation equation for a sampling data stream to be supplied since time t+1 wiil be described below. 

[0082] That is, the sampling data stream supplied since the time t+1 is expressed as a set of x(t+1 ), x(t-i-2), x{t+3) 

30 x(t+N-1), x(t+N). Because data x(t) is deleted from data stream supplied since time t while new data is added, trans- 
formation of the real part is expressed in a following equation and that expressed equation is developed as follows. 
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[0083] Then, if the oldest data x(t) is deleted and the newest data x(t+N) is incorporated, Y^Ck, t+1 ) can be expressed 

with Xr(k. t) and likewise, Yj(k. t+1) can be also expressed with Xj(k, t). 

[0084] As a result, the equation (6) and equation (7) are expressed as follows. 
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[0085] Next, the transformation equation for the imaginary part is also developed in the sanne way. 

[0086] That is, when new data x(t+N) is supplied, the transformation of the imaginary part is also developed as shown 

in a following equation. 
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[0087] Therefore, the equations (5), (10) can be expressed based on the relation of the equations (8), (9). 
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[0088] In this way, the discrete Fourier transformation at time t+1 can be executed according to the equations (11), 
(12) by using the result of the discrete Fourier arithmetic operation at time t recursively and further using data value 
which is a difference between newly fetched sample value x(t+N) and a sample value x(t) to be deleted. 
[0089] It means that a result of Fourier transfomnation on the newest N-point data containing new data obtained by 
successive sampling at the sampling time interval shown in Fig. 2 can be obtained. 

[0090] Here, if it is regarded that the Fourier transfomnation executed in the above-described manner is a linearly 
invariable system in which continuous sampling data is supplied to a system having a certain coefficient and converted, 
this conversion processing can be regarded lo be filtering processing. 

[0091] Then, if 1/(square root of N) which gives an amplitude is assumed to be A and the real part Xr(k,t) and the 
imaginary part Xi(k,t) of Fourier coefficient with respect to data stream supplied since time t are assumed to be y^^ (t) 
and yi,^ (t) respectively, the aforementioned equations (11), (12) are treated as a transformation equation of output y 
to input X, so that the transfomnation equation of the Fourier coefficient of data supplied since time t+1 can be expressed 
as follows. 
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y,^ + 1) = {y,^ (? ) + A[x{t + jV ) - x{t )]}x COS [2^1 + y.^ {t ) sin 



N 



(13) 



ya(f + !) = >',.* (0 COS 



2^^ 



- {y ... ) + ^ + JS^ ) - )]}x sin 



2^ 
N 



•••(14) 
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(where A = —j=^ ) 

[0092] Further, if these equations (13), (14) are subjected to Z transformation, a result of the transformation can be 
expressed as shown in following equations (15), (16). 



20 



y,,(z)= ^A[x{zy -X{z)]xcos [2^1 +I^,,(2)r-'sin 



2* 

N 



(15) 



n*(^) = i^.*U>"' COS 
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-f,,,(z)z-' +A[x{zy -A'(z)|sin 



(16) 



[0093] If a transfer function of output Y(2) to input X(2) is assumed to be H(z), H(z) is given in a following equation. 

H(z) = Y{z)/X{z) 

[0094] Then, if the equations (15), (16) are treated as simultaneous equations and it is Intended to obtain each 
transfer function of the real part and imaginary part, they can be obtained as shown in following equations (17), (18). 



COS 


2*1 

N 




1-2 cos 




N 


z-'^z-^ 



- (17) 
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HXz) = a{i-z-'\ 
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-sin 


N 




1-2 cos 


N 


z~ 





(18) 



[0095] Like this , a transfer function indicating the recursive DFT is expressed in a following equation (1 9). 
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[0096] In this equation, t-Z'^ indicates deletion of the oldest data upon updating the Input data and a portion ex- 
10 pressed with a fraction is so constructed that its numerator and its denominator are IIR filters (cyclic filter) having 
- primary and secondary degrees respectively. 

[0097] Although the DFT is a FIR filter for carrying out processing to a finite point number N. if it is achieved with a 
MR filter as shown in the aforementioned equation, it can be constructed with a simpler hardware than the FIR filter 
[0098] Fig.3 shows the structure for carrying out the discrete Fourier transfonnation achieved with the IIR filter 
75 [0099] In the same Figure, data sampled at a predetemnined sampling ciock, for example, data x(t+N) sampled at 
time t+N is supplied to a Fourier operating portion 5. Here, the supplied sampling data Is computed by the characteristic 
of a transfer function indicated as H(z). so that X^(K t+1 ) and jXj (k, t+1 ) are supplied as real part output and imaginary 
part output of the Fourier coefficient. 

[0100] The Fourier operating portion 5 having the transfer function H(z) indicated by the equation (19) carries out 
20 IIR filtering operation in succession for every supplied sampling data and supplies results of the operation of the DFT 
transformation corresponding to the real part and the imaginary part. 

[0101] That is, if data sampled at a certain sampling time t+l is supplied to the Fourier operating portion 5, complex 
Fourier coefficient X^ik, t+1)-jXj(k,t+1) corresponding to N sampling data supplied since the sampling time t+1 can be 
obtained. 

25 [0102] As described above, the Fourier operating portion 5 can carry out Fourier transfonnation processing on con- 
tinuously supplied data stream continuously at the sampling time interval. Then, a result of arithmetic operation can 
be obtained by changing analysis resolution in the Fourier transfonnation, so that high-resolution frequency analysis 
can be carried out in time direction by changing the length N of data stream supplied, for example, increasing the value 
of N. 

30 [0103] As for the change of the analysis resolution, in case of adaptive signal processing to be carried out using a 
result of analysis by, for example, the DFT operation or adaptive control for control system, which is an object, the 
discrete Fourier transfonnation needs to be analyzed to arbitrary frequency band at an arbitrary resolution (frequency 
interval). 

[0104] Fig.4 shows the relation of frequency, which should be obtained by an object discrete Fourier analysis. 

35 [0105] The same Figure shows recursive discrete Fourier transformation in which the minimum frequency of a fre- 
quency band to be analyzed is f 1 while the maximum frequency thereof is f2 and N-point resolution is provided between 
f 1 and f2- Then, a method of carrying out the DFT under such indicated resolution will be described. 
[0106] Setting of the resolution is carried out by the basic frequency setting portion 2 represented on the recursive 
discrete Fourier transfonnation device shown in Fig. 1 and based on the resolution infonnation set up here, the Fourier 

40 transformation is carried out such that while the minimum frequency of arbitrary frequency band is f 1 and the maximum 
frequency is f2, a desired resolution (frequency interval), indicated with any of values which divides an interval between 
f1 and f2 to N, is N point. 

[0107] A method of carrying out the discrete Fourier transfonnation operation with a desired resolution in a desired 
arbitrary frequency band at every sampling time interval will be described in detail. 
45 [0108] In the equations (11) and (12), or (19) used for carrying out the recursive discrete Fourier transformation, cos 
[2 nk/N] and sin[27ck/N] expressed with circular functions indicate basic functions relating to analysis point nunnber N 
which should be obtained by operation. 

[0109] Therefore, the values of f1 , f2 and fs are set up so that k/N relating to the basic frequency indicated in the 
circular functions coincides with a desired frequency band. 
50 [0110] Here^ it is assumed that basic frequency a of a basic function corresponding to the minimum frequency f1 is 
a = fl/fs and basic frequency p of a basic function corresponding to the maximum frequency f2 is p = f2/fs. 
[0111] Then, when the resolution is N points, the basic frequency in the frequency band of f1 to f2 is as follows. 
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(where k = 0 - N-1) 

Therefore, by replacing the basic frequency k/N in the equations (11) and (12) or (19) with a following equation (20), 
it is possible to obtain recursive Fourier transformation equation adaptive for transformation with a desired arbitrary 
resolution in an arbitrary frequency band desired to be analyzed. 



1 r(£Wi). 1 

fs\ N-1 j 



(20) 



(where 0<,k<N- ^) 

[01 1 2] Equations for carrying out transformation at an arbitrary frequency resolution, wliich corresponds to the equa- 
tions (11), (12) are expressed as equations (21), (22). 



■(/2-/l)fe 



X COS i 2 



71 



fs 



N-1 



+ /1 



{f2-fl)k 



N-1 



+ /1 



•• (21) 



30 



35 



40 



45 



{f2-n)k 



N-1 



X COS J 2 — 



(22) 



[0113] An equation for carrying out transfonnation at an arbitrary frequency resolution, which corresponds to the 
equation (19), is expressed as an equation (23). 



50 



/f(z) = ^(l-Z-< 



cos[2j57]- sin 




l-2cos[2;57 





(23) 



55 



where A is a positive constant for providing [x(t+N) - x(t)] with an annplitude value, and 
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[0114] Fig. 5 shows the structure for carrying out the Fourier transformation at an arbitrary frequency resolution. 
[0115] In the same Figure, sampling data sampled at every equal time interval is supplied to a Fourier operating 
portion 6, which is a linearly invariable system, in succession, so that for example, data x(t+N) is supplied at time t+N. 

10 [0116] At this time, a basic frequency setting portion 7 sets up a lower limit frequency f1 and a upper limit frequency 
f2 of an arbitrary frequency band to be analyzed and then a desired resolution N. Then, a basic frequency to be trans- 
formed is determined according to the aforementioned equation (20) and supplied to the Fourier operating portion 6, 
which is a linearly invariable system. The Fourier operating portion 6 carries out IIR filter processing having a transfer 
function indicated by the aforementioned equation (1 9) according to the basic frequency. 

^5 [0117] The Fourier operating portion 6 supplies Kj{k, t+1) as the real part of the Fourier coefficient to N sampling 
data supplied since time t+1 and Xj(k, t+1) as the imaginary part thereof. Consequently, the DFT arithmetic operation 
can be carried out in a desired arbitrary frequency band at a desired resolution at every sampling time interval. 
[0118] As described above, recursive discrete Fourier transformation operation, in which the resolution is set up 
arbitrarily for analysis, can be executed. This Fourier transformation processing can be candied out based on the equa- 

20 tions (11) and (12), or equation (19) for every sampling value supplied newly. By changing the value of k to 0, 1 , 2, .... 
N-1 in that operation, a result of transformation processing under each basic frequency set corresponding to each k 
is provided. 

[0119] Because in the transformation equation for conducting recursive discrete Fourier transformation operation, 
the basic frequency set by a predetermined k value is capable of obtaining a result of transformation processing inde- 
25 pendently of analysis on the basic frequency set by other k value, a result of Fourier transformation operation on a 
specific basic frequency relating to ariDitrary k value may be obtained independently. 

[0120] Then, because when attention is paid to a specific frequency point or k is assumed to be a specific value, the 
value of circular functions in the equations (11) and (12). or equation (19) is a constant. Thus, by using the aforemen- 
tioned definition equation (1) for the discrete Fourier transfomnation, which has been employed since before, a result 
30 of the Fourier operation relating to that basic function can be obtained more easily than the Fourier transformation 
relating to a specific basic function. 

[0121] Because the arithmetic operation for the signal component relating to the specific basic function is more 
difficult in case of the FFT which uses butterfly operation, the recursive discrete Fourier transfomnation method shown 
here is expected to be widely applied as a signal analysis method, which will be employed in more diversified ways in 



[0122] The recursive discrete Fourier transformation method has been described above. Because this transformation 
method is capable of selecting a basic frequency preferable for analysis, the number of points to be handled is not 
restricted to power of 2 unlike in case of the FFT, so that any positive integer can be adopted. 
[0123] Although in the above description, it is assumed that the sampling data to be supplied is real number data, 
40 even if the sampling data is complex data like shown below, it can be expanded. 

[0124] Specifically, in case where the sampling data is complex data, an equation corresponding to the equations 
(3), (4) is as follows. 
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^,(A',/ + l)=|z,(fc,0+^kO+A^)-x,(t)]|cos[2^ 



|A:,M+^[x,{/+A0-x,(r)]|sin[2^] - (25) 



(where, >4= jN,0<k<N- 1) 
[0125] Next, inverse discrete Fourier transformation, which is a second ennbodinnent achieved by applying the above 
discrete Fourier transformation technology, will be described below. 

[0126] Fig. 6 shows an embodiment of the recursive inverse discrete Fourier transformation device, to which the 
15 recursive Inverse discrete Fourier transformation method is applied. This embodiment will be described with reference 
to the same Figure. It is assunried that the sarnpling data to be supplied is complex data. 

[0127], This recursive inverse discrete Fourier transfomnafibn device compnses.a data updating portion 11 , which is 
supplied with data sampled at every predetermined time interval, for storirig^iaiDpTfed new N-data (N: positive integer) 
temporarily, a basic frequency setting portion 12 for selling the basic frequency for carrying oul the Inverse discrete 
20 Fourier transformation, a recursive inverse discrete Fourier transformation (ID FT) operating portion 13 and a memory 
portion 14 tor storing computed data temporarily. 

[0128] An operation of the recursive inverse discrete Fourier transformation device having such a structure will be 
described. 

[0129] First, supplied complex data is sampled by a sampling circuit (not shown) at every predetemnined time interval, 
25 so that discrete data of the real part and imaginary part, sampled and quantitized, are supplied to the data updating 
portion 11 . 

[0130] The sampling circuit samples data values supplied at times t, t+1 , t-t-2, t+3 t+N-1 , t+N (N: natural number) 

each having an equal interval so as to generate data stream x^{t) + jXj, Xr(t+1) + jXj(t+1), x^(t+2) + jXj(t-i-2), Xr(t+3) + jXj 
(t+3), Xr(t+N-1 ) + jXi(t+N-1), Xr(t+N) +jXj(t-t-N) with data value supplied at each time as sampling value corresponding 
30 to each time. 

[0131] In this manner, the sampling data Xr(t) + jXi(t) sampled at time t, Xr(t+1) + jXj(t+1) at time t+1 Xr(t+N-1) + 

jXi(t+N-1 ) at time t+N-1 Xr(t+N-1 ) + jXj(t+N-1 ) at time t+N-1 , Xr(t+N) + jXj(t+N) at time t+N .... are supplied to the data 

updating portion 11. The data updating portion 11 updates complex data of new N data (N: positive integer) of the 
supplied data and stores them therein temporarily. 

35 [01 32] The N data stored temporarily is supplied to the recursive ID FT operating portion 1 3 and a last FFT operation 
result stored in the memory 14 temporarily is supplied to the recursive IDFT operating portion 13 as recursive data. 
Then, the recursive IDFT operating portion 13 carries out recursive discrete inverse Fourier transformation according 
to frequency resolution infonnation set up by the basic frequency setting portion 12 in a manner described later and 
outputs a result of Its arithmetic operation. 

40 [0133] Then, the operation result subjected to inverse discrete Fourier transformation by operation by the recursive 
IDFT operating portion 13 is supplied to the memory portion 14 and stored therein temporarily. The temporarily stored 
operation result is supplied to the recursive IDFT operating portion 1 3. Here, the IDFT operation for next supplied data 
is carried out, so that the result of the IDFT operation is used recursively. 
[0134] Next, the IDFT operation, which is carried out recursively, will be described in detail. 

45 [0135] Inverse Fourier transformation X(k,t) at arbitrary time t is defined according to a following equation (26). 
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Sin 



2{n-t) — 



10 



4 ;J-7= y.r,(n)sin 



2{n-t) 



= ^x,{n)cos 2{n-t) 



»5 



for k= 0,1. ..... N-1 

20 [0136] Data stream x(n) supplied here is complex number which is a number expressed like x(n) = x^(n) + jXj(n). 

[0137] The recursive inverse discrete Fourier transformation is carried out for N complex data supplied in this manner 
and its arithmetic operation is carried out through the same processing as the aforementioned recursive DFT transfor- 
mation method. 

[0138] In this operation, complex Fourier coefficient Xf(k, t+1) -i- jXj(k, t+l) at sampling time tn 1 can be obtained 
25 according to following equations (27), (28) by using a difference between newly complex dataXr(t+N) + jXj(t+N) sampled 
and fetched at sampling time t+1 and the oldest complex data Xr(t) + jX((t) of N-data stream and complex Fourier 
coefficient X^{k, t) + jXj(k; t) already obtained at sampling time t. 
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35 



X,{k,0.^M.N)-xMYin^2f 



(27) 



40 



X,{k,t +l)= |;^r,(/:,r) + -i-[x,(/ + Ar)-x,(/)]j 



cos 



2^ 

N 



(28) 



so (where, B^jN,0<k<N- 1) 

[0139] The recursive IDFTtransfomnation is carried out in the above way. Like the Fourier transfonnation, the inverse 
Fourier transfomnation is also capable of transforming the newest N-point data stream containing new data, which is 
sampled successively at every sampling time interval shown in Fig. 2. 

[0140] Meanwhile, in case where the sampling data is complex data, an equation corresponding to the equations 
55 (21), (22) is as follows. 
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;...(.,r.l) = {x.(.,0^^k(..^)-...{0]}xcosj2^[%i^^ 



(30) 



[01411 corresponding to the equations (29), (30). an equation for carrying out inverse transfomnation at an arbitrary 
frequency resolution is expressed in the form of equations (31), (32). 
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[0142] It should be understood that many modifications and adaptations of the invention v^ill become apparent to 
ihose Skilled in the art and il is intended lo encompass such obvious modifications and changes .n the scope of the 
claims appended hereto. 

Claims 

1. ArecursivediscreteFouriertransformationmethodwhereindatavaluesx(t),x(t4.1).(t.2),x(t.3), xaj-l^^ 

sampled at times t, t+1 , t^2. t+3 t+N-1 , t+N (N is a positive integer which .s 1 or more) each having an equal 

interval are supplied and with such N data values supplied since time t as a data stream, afrequency cornponen , 
wwch is degree k (k is 0 or a positive integer smaller than N) obtained by carrying out complex ^ouner ansio. 
mation on the data stream, is obtained such that a real part X^k, t) thereof and an imaginary part X,(k, t) thereof 
are complex Fourier coefficients, the method comprising: 



22 



30CID- <EP 1168192A2J_> 



EPfl ^ 192 A2 



10 



a first step of storing the data stream x(t), x(t+1), x(t+2), x(t+3), ... x(t+N-1) supplied since tinne t at time t+N- 
1 into a first merfiory means (1) temporarily; 

a second step of obtaining the complex Fourier coefficients Xf(k, t) and Xj(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients X^{k, t) and Xj(k, t) obtained at the second step into the 
second memory means (4) temporarily; and 

a fourth step of by using data value x(t+N) supplied at time t-i-N, data value x(t) stored in the first memory 
means (1 ) temporarily and the complex Fourier coefficients X^{K t), Xj(k,t) stored in the second memory means 
(4) temporarily, obtaining complex Fourier coefficients Xr(k, t+1) and Xi(k, t+1) to the data stream supplied 
since time t+1 with respect to a positive constant value A for giving an amplitude value corresponding to a 
difference between the x(t+N) and the x(t) according to a following equation: 



15 



X,{k,t^l) = S^X^ {k,t ) + j-[x{t + N)- x{t )]| 



X COS 
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X,{k,t +l) = Xi{k,t)cos 



2fYi^xAk,t).^[x{t.N)-4t)] 







.sin 
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2. A recursive discrete Fourier transformation method according to claim 1 wherein inverse discrete Fourier trans- 
formation data subjected to inverse discrete Fourier transformation using a positive constant a is supplied so as 
to obtain a complex Fourier transformation coefficient to the supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out using the positive constant A corresponding to the pos- 
itive constant a. 

3. A recursive discrete Fourier transformation method according to claim 1 wherein k indicating a degree of the 
complex Fourier coefficient is set to a desired value and Fourier transfonmation to the set value k is repeated at 
every time having an equal interval. 
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A recursive discrete Fouriertransfomnation method wherein with sampling frequency as fs, data values x(t), x(t+1), 
x{t-i-2), x(t+3), .... x{t+N-1), x(t+N) sampled at times 1. 1+1 , t+2. t+3, t+N-1 , t+N (N is a positive integer which is 
1 or more) each having an equal interval are supplied and with such N data values supplied since time t as a data 
stream, a frequency interval given with a minimum frequency f1 and a maximum frequency f2 to the data stream 
is regarded as a measuring frequency interval, a frequency interval obtained by dividing that measuring frequency 
interval by the N is assumed to be an analysis frequency interval and a result of the frequency analysis obtained 
by carrying out complex Fourier transformation at every analysis frequency interval is obtained as a frequency 
component which is k (k is 0 or a positive integer smallerthan N) times the frequency interval, such that a real part 
Xr(k, t) thereof and an imaginary part Xj(k, t) thereof are complex Fourier coefficients, the method comprising: 

a first step of storing the data stream x(t), x(t+1 ), x(t+2), x(t+3) x{t+N-1 ) supplied since time t at time t+N- 

1 into a first memory means (1) temporarily; 

a second step of obtaining the complex Fourier coefficients Xf(k, I) and Xj(k, I) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients Xr(k, t) and Xj(k, t) obtained at the second step into the 
second memory means (4) temporarily; and 

a fourth step of by using data value x(t+N) supplied at time t+N, data value x(t) stored in the first memory 
means (1 ) temporarily and the complex Fourier coefficients X^{k, t), Xj(k,t) stored in the second memory means 
(4) temporarily, obtaining complex Fourier coefficients X^{k, t+1) and Xj(k, t+1) within the frequency interval 
with the minimum frequency f 1 and the maximum frequency f2 to the data stream supplied since time t+1 with 
respect to a positive constant A for giving an amplitude corresponding to a difference between the x(t+N) and 
the x(t) according to following equations: 
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X sin J 2 
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5 A recursive discrete Fourier transformatior, method according to claim 4 wherein inverse discrete Fourier trans- 
formation data subjected to inverse discrete Fourier transformation using a positive constant a ,s supplied so as 
to obtain a complex Fourier transfomiation coefficient to the supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out using the positive constant A corresponding to the pos- 
itive constant a. 

6 A recursive discrete Fourier transformation method according to claim 4 wherein k indicating a degree of the 

■ complex Fourier coefficient Is set to a desired value and Fourier transformation to the set value k is repeated at 
every time having an equal Interval. 

7 A recursive discrete Fourier transformation method wherein data values x(t),x(t-i-1),x(t-i-2),x(t-i-3), .. x(t+N-1). x 

■ (t+N) sampled at times t, t+1 , t+2, t+3 t+N-1 , t-.N (N is a positive integer which is 1 or more) each having an 

equal inteLl are supplied and with such N data values supplied since time t as a data stream , a ''^^^l^^^^^^ 
ponent, which Is degree k (k is 0 or a positive integer smaller than N) obtained by carrying out complex Foune 
transfomiation on the data stream, is obtained such that a real part X,(k, t) thereof and an imaginary part X^K t) 
thereof are complex Fourier coefficients, the method comprising: 

a first step of storing the data stream x(t), x(t+1), x(t^2), x(t+3) x(t-^N-1) supplied since time t at time t+N- 

1 into a first memory means (1 ) temporarily; ^ ^ ^ , . ♦^^^w 

a second step of obtaining the complex Fourier coefficients X,(k, t) and Xi(k, t) of the data stream stored 

temDorarily at the first step; ^ . . „ 

a third step of storing the complex Fourier coefficients X,(k, t) and Xi(k, t) obtained at the second step into the 

45 second memory means (4) tennporarily; and ^,^«ono 

a fourth step of by using the complex Fourier coefficients X,(k. t) - jX^k. t) stored in the second memory means 
(4).temporarily, obtaining complex Fourier coefficients X,(k, t+1 ) - M ) to the data stream supplied since 
time t+l based on a transfer function expressed in a following equation. 
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Where A is a positive constant for providing [x(t+N) - x(t)] with an amplitude. 
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8. A recursive discrete Fourier transformation method according to claim 7 wherein the positive constant value A for 
providing with an amplitude corresponding to a difference between the x(t+N) and the x(t) is capable of being set 
selectively with an inverse number of square root of 1 , N or 1/N. 

5 9. A recursive discrete Fourier transformation method according to claim 7 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the Fourier transformation to the set value k is repeated 
at every tinne having an equal interval. 

10. A recursive discrete Fourier transformation method wherein with sampling frequency as fs, data values x(t), x (t+1), 
10 x(t+2). x(t4-3), .... x(t+N-1), x(t+N) sampled at times t, t+1 , t+2. t-»-3, .... t+N-l , t+N (N is a positive integer which is 

1 or more) each having an equal interval are supplied and with such N data values supplied since time t as a data 
stream, a frequency interval given with a minimum frequency f1 and a maximum frequency f2 to the data stream 
is regarded as a measuring frequency interval, a frequency interval obtained by dividing that measuring frequency 
interval by the N is assumed to be an analysis frequency interval and a result of the frequency analysis provided 
15 by carrying out complex Fourier transfomnation at every analysis frequency interval is obtained as a frequency 

component which is k (kis 0 or a positive integer smaller than N) times the frequency interval, such that a real part 
Xr(k, t) thereof and an imaginary part Xi(k, t) thereof are complex Fourier coefficients, the method comprising: 

a rirsl step of storing the data stream x(t), x(l+1). x(l-i-2), x{l+3), x(l+N-1) supplied since time I al lime l+N- 
20 1 into a first memory means (1) temporarily; 

a second step of obtaining the complex Fourier coefficients K^(k, t+l) - jXj(k, t-j-1) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients X^(k, t+1) - jXj(k, t+1) obtained at the second step into 
the second memory moans (4) temporarily; and 
25 a fourth step of by using the complex Fourier coefficient Xf(k, t) - jXj(k, t) stored in the second memory means 

(4) temporarily, obtaining complex Fourier coefficients X^Ck, t+1 ) - jX,(k, t+1 ) at the frequency interval given by 
the minimum frequency f1 and the maximum frequency f2 to the data stream supplied since time t+1 based 
on a transfer function expressed in a following equation: 
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COS [Tjz p]- j sin [Ztt p]- z' 



l-2cos[2;r p]z~^ + z" 
where A is a positive constant for providing [x(t+N) - x(t)] with an amplitude value, and 



11. A recursive discrete Fourier transfomnation method according to claim 10 wherein the positive constant value A 
45 for providing with an amplitude corresponding to a difference between the x(t+N) and the x(t) is capable of being 

set selectively with an inverse number of square root of 1 , N or 1/N. 

12. A recursive discrete Fourier transfomnation method according to claim 10 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the Fourier transformation to the set value k is repeated 

50 at every time having an equal interval. 

13. A recursive discrete Fourier transformation method for obtaining complex Fourier coefficients, wherein data values 
\{\) + jXj(t), x,(t+1 ) + jXj{t+1 ), Xr(t+2) + jXj(t+2), x,(t+3) + jXj(t+3). x,(t+N-1 ) + jXi(t+N-1 ), x^Ct+N) + jXj(t+N) sampled 
at times t, t+1 , t+2, t+3 , t+N-1, t+N (N is a positive integer which is 1 or more) each having an equal interval 

55 are supplied and with such N data values supplied since time t as a data stream, a frequency component, which 

is degree k (k is 0 or a positive integer smaller than N) obtained by carrying out complex Fourier transformation 
on the data stream, is obtained such that a real part Xf(k, t) thereof and an imaginary part Xj(k, t) thereof are 
complex Fourier coefficients, the method comprising: 
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a first step of storing the data stream x,(t) + jXj(t), x,(X+1) + jXj(t+1)/xr(t+2) + jXj(t+2), x,(t+3) + jXj(t+3), x, 
(t+N-1) + jXj(t+N-1) supplied since time t at time t+N-1 into a first memory means (1) temporarily; 
a second step of obtaining the complex Fourier coefficients Xr(k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients X,(k, t) and Xj(k, t) obtained at the second step into the 
second memory means (4) temporarily; and 

a fourth step of by using the data value Xr(t+N) + jXi(t+N) supplied at time t+N, the data value x,(t) + jXi(t) stored 
in the first memory means (1 ) temporarily and the complex Fourier coefficient Xr(k, t) and Xj(k, t) stored in the 
second memory means (4) temporarily and used recursively, obtaining complex Fourier coefficients X^(k, t+1 ) 
and Xj(k, t+t) to the data stream supplied since time t+1 with respect to a positive constant A for giving an 
amplitude value corresponding to the difference between the \(t+N) and the x^O) according to a following 
equation: 
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^(A:,^+l) = |z,(/:,0 + ^k(^ + ^)-x,(^)]| 
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14. A recursive discrete Fourier transformation method according to claim 13 wherein by using data value y^{t) + jyj(t) 
supplied at timet, data value y^lt-i-N) +jyj(t+N) supplied at time t+N and complex inverse discrete Fourier coefficients 
Yr(k, t) and Yi(k, t) with the real part and the imaginary part obtained with respect to N data values y^(\) + jyj(t), 
(t+1) +jyi(t+1), .... yr(t+N-1)+jyi(t+N-1) supplied from timet to time t+N-1, 

a real part t+1) and an imaginary part Yj(k, t+1) of each complex inverse discrete Fourier coefficient of 
N data values supplied since time t+1 are obtained with respect to a positive constant value B for giving an 
amplitude corresponding to a difference between the yr(t+N) and the y^lt) as inverse discrete Fourier transfor- 
mation data, according to following equations, 
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and afterthe obtained inverse discrete Fourier transformation data is supplied, discrete Fourier transfomnation 
on the supplied inverse discrete Fourier transfomnation data is carried out. 

the recursive discrete Fourier transformation conducting discrete Fourier transformation by using the positive 
constant value A corresponding to the positive constant value B. 

15. A recursive discrete Fourier transfomnation method according to claim 13 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the Fourier transfonmation to the set value k is repeated 
at every time having an equal interval. 



16. 



A recursive discrete Fourier transformation method for obtaining complex Fourier coefficients, wherein with sam- 
pling frequency as fs. data values Xr(t) + jXj(t), Xr(t+1) + jXj(t+1), Xr(t+2) + jXj(t+2), Xr(t+3) + jXj(t+3), Xr(t+N-1) + 
jXj{t+N-1), Xr(t+N) + jXj(t+N) sampled at times t, t+1 , t+2, t+3, t+N-1 , t+N (N is a positive integer which is 1 or 
more) each having an equal interval are supplied and with such N complex data values supplied since time t as a 
data stream, a frequency interval given with a minimum frequency f1 and a maximum frequency f2 to the data 
stream is regarded as a specified frequency interval, a frequency interval obtained by dividing the specified fre- 
quency interval by the N is assumed to be a minimum frequency interval and a result of the Fourier transformation 
provided by carrying out complex Fourier transformation at every minimum frequency inten/al is obtained as a 
frequency component which is k (k is 0 or a positive integer smaller than N) times the minimum frequency interval, 
such that a real part X^(k, t) thereof and an imaginary part Xj(k, t) thereof are complex Fourier coefficients, the 
method comprising: 

a first step of storing the data stream + jXi(t), x,(X+^) + jXj(t-i-1), Xr(t+2) + jXj(t+2). Xr(t+3) + jXi(t+3), x^ 
(t+N-1) + jXi(t+N-1), Xr(t-i-N) + jXj(t+N) supplied since time t at time t+N-1 into a first memory means (4) tem- 
porarily; 

a second step of obtaining the complex Fourier coefficients X^Ck, t) and Xj(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients Xr(k, t) and Xj(k, t) obtained at the second step into the 
second memory means (4) temporarily; and 

a fourth step of by using the data value Xr(t+N) + jXi(t4-N) supplied at time t+N, the data value Xr(t) + jXj(t) stored 
temporarily in the first memory means (1 ) and complex Fourier coefficients X^(k, t) and Xj(k, t) stored temporarily 
in the secondmemorymeans (4) and used recursively, obtaining complex Fourier coefficient X^lk, t+1) and Xj 
(k, t+1 ) in the specified frequency interval given by the minimum frequency f 1 and the maximum frequency f2 
to the data stream supplied since time t+1 with respect to a positive constant A for giving an amplitude value 
corresponding to a difference between the x(t+N) and the x(t) according to a following equation: 
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1 7. A recursive discrete Fourier transformation method according to claim 1 6 wherein by using data value y^O) + jyj(t) 
supplied at timet, data value y^U+N) + jyi(t+N) supplied at time t+N and complex inverse discrete Fourier coefficients 
^5 Yr(k, t) and Yj(k, t) with the real part and the imaginary part obtained with respect to N data values y,{\) + jyj(t), y 
(t+1) + jyi(t+1) y^(t-HN-l) + jyj(t+N-1) supplied from time t to time t+N-1 , 

a real part Y,(k, t+1) and an imaginary part Yi(k, t+1) of each complex inverse discrete Fourier coefficient of 
N data values supplied since lime 1+1 are obtained with respect to a positive constant value B for giving an 
amplitude corresponding to a difference between the y^Ct+N) and the y,(t) as inverse discrete Fourier transfor- 
mation data, according to following equations, 
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and afterthe obtained inverse discrete Fourier transformation data is supplied, discrete Fourier transformation 
on the supplied inverse discrete Fourier transfonnation data is carried out, 

the recursive discrete Fourier transformation conducting discrete Fourier transformation by using the positive 
constant value A corresponding to the positive constant value B. 

18. A recursive discrete Fourier transformation method according to claim 16 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the Fourier transformation to the set value k is repeated 
at every time having an equal interval. 



19 



A recursive inverse discrete Fourier transformation method for obtaining complex inverse Fourier coefficients 
wherein complex data values x^^) + jx,(t), x,(t+1 ) + jx,(t+1 ), x,(t+2) + jXi(t+2). x,(t+3) + jx,(t+3), x,(t+N-1 ) 4-jx,(t+N- 
1), Xr(t+N) + jXi(t+N) sampled at times 1. 1+1 , t+2, t+3, t+N-1 . t+N (N is a positive integer which is 1 or more) 
each having an equal interval are supplied and with such N data values supplied since time t as a data stream, 
an inverse complex Fourier transfonnation component, which is degree k (kis 0 or a positive integer smaller than 
N) obtained by carrying out inverse complex Fourier transfonnation on the data stream, is obtained such that a 
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real part \{k, X) thereof and an imaginary part Xj(k. t) thereof are inverse complex Fourier coefficients, the method 
comprising: 

a first step of storing the data stream x,(t)+jXi(t), x^(t+1 ) + jXi(t+1 ), x^(t+2) + jXi(t+2), x^(t+3) + jXi(t+3). x^t+N- 
1) + jXj(t+N-1) supplied since time t at timet+N-1 into a first memory means (11) temporarily; 
a second step of obtaining the inverse complex Fourier coefficients X^Ck, t) and Xi(k. t) of the data stream 
stored temporarily at the first step; 

a third step of storing the inverse complex Fourier coefficients Xj(k, X) and Xj(k. t) obtained at the second step 
into the second memory means (14) temporarily; and 

a fourth step of by using the data value Xr(t+N) + iXj(t+N) supplied at time t+N, the data value x^(X) + jXj(t) stored 
in the first memory means (11) temporarily and the inverse complex Fourier coefficient X/k, t) and Xi(k, t) 
stored in the second memory means (1 4) temporarily, obtaining inverse complex Fourier coefficients X^(k, t+1 ) 
and Xj(k, t+1) to the data stream supplied since time t+1 with respect to a positive constant B for giving an 
amplitude value corresponding to the difference between the Xr(t+N) and the Xr(t) according to a following 
equation: 
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20. A recursive inverse discrete Fourier transformation method according to claim 1 9 wherein k indicating the degree 
40 of the complex Fourier coefficient is set to a desired value and the inverse Fourier transformation to the set value 

k is repeated at every time having an equal interval. 

21. A recursive inverse discrete Fourier transformation method for obtaining complex inverse Fourier coefficients 
wherein with sampling frequency as fs, data values Xr(t) + jXj(t), Xr(t+1) + jXi(t+1), Xr(t+2) + jXi(t+2), Xr(t+3) + jXj 

45 (t+3), Xr(t+N-1) + jXi(t+N-1), Xr(t+N) + jXj(t+N) Sampled at times t, t+1, t+2, t+3, .... t+N-1, t+N (N is a positive 

integer which is 1 or more) each having an equal interval are supplied and with such N data values supplied since 
lime I as a data stream, a frequency interval given with a minimum frequency f1 and a maximum frequency f2 to 
the data stream is regarded as a specified frequency interval, a frequency interval obtained by dividing that specified 
frequency interval by the N is assumed to be a minimum frequency interval and a result of the inverse Fourier 

50 transformation provided by carrying out inverse complex Fourier transfonnation at every minimum frequency in- 

terval is obtained as a frequency component which is k (k is 0 or a positive integer smaller than N) times the 
minimum frequency interval, in the form of a real part X^(k. t) and an imaginary part Xj{k, t). the method comprising: 

a first step of storing the data stream Xr(t) +jx,(t), Xr(t+1) + jx,(t+1), Xr(t+2) + jx,(t+2). Xr(t+3) + jx,(t+3) x^(t+N- 
55 1 ) + jXj(t+N-1 ), Xr(t+N) + jXj(t+N) supplied since time t at time t+N-1 into a first memory means (1 1 ) temporarily; 

a second step of obtaining the inverse Fourier coefficients X^{k, t) and Xj(k, t) of the data stream stored tem- 
porarily at the first step; 

a third step of storing the inverse Fourier coefficients Xr(k, t) and Xj(k, t) obtained at the second step into the 
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second memory means (14) temporarily; and 

a fourth step of by using the data value Xr(t+N) + jXj(t+N) supplied at time t+N, the data value Xr(t) + jXj(t) stored 
temporarily in the first memory means (11 ) and preceding inverse Fourier coefficients X,{k, t) and Xi(k, t) stored 
temporarily in the second memory means (14), obtaining inverse Fourier coefficient X^{k, t+1) and Xj(k, t+1) 
in the specified frequency interval given by the minimum frequency f1 and the maximum frequency f2 to the 
data stream supplied since time t+1 with respect to a positive constant B for "giving an amplitude value corre- 
sponding to a difference between the x(t-i-N) and the x(t) according to following equations: 
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X, {k,i + 1) = |;ir, + + at) - {t )]| x cos |2^ 
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22. A recursive inverse discrete Fourier transfonnation method according to claim 21 wherein k indicating the degree 
of the complex Fourier coefficient is set to a desired value and the inverse Fourier transformation to the set value 
k is repeated at every tinne having an equal interval. 



35 



40 



45 



50 



55 



aa,. 

OOCID: <EP 1168192A2_L> 



H 
D 
O 



+ VII 
+ VII 



O 

I— ^ 

tin 





EQUENCY 
PORTION 


CO— ^ 


TION 




> 


r 






RECURSIVE DFT 
OPERATING POR 




[ORY 











MEM 




< tu 

P3 CO 












X 



5iXX;iO: <EP 1169192A2J_> 



31 



EP 1 168 19fjA2Wai 




»DCX;iD: <EP 1l6ei92A2J_> 



32 



EP 



d 
I— I 




Q 
O 

z 

DO 



+ 

X 



SI 
I 



O 
O 



CM 

I 

+ 



to 



CO 

o 

a 



2 

I 

< 



SDOCID: <EP 1 168192A2_L> 



33 



EP 1 168 igi^A24^t 



FIG.4 



LEVEL 



r 



RESOLUTION LEVEL 

N POINT 
A 



fl 



f2 



34 



EP ;|}.^ 1^2 A2 



I 



I 



H 

D 

H ^ 
ID ^ 
O X 



< 




< 

o 
< 



ONI. 




w 
O 


X 


FOURIER 
PORTION 





Q 
U 

< 
CO 



X 



35 

SDCX;lD: <EP 1168192A2_L>- 



EP 1 168 ^99)A^Jk*i 



H 
O 



T VII 
VII 

X 




iDOCID: <EP 1168192A2J_> 



36 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ BVlAGE CUT OFF AT TOP, BOTTOM OR SIDES 
ID FADED TEXT OR DRAWING 

0 BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: \ 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, pleajse do not report these problems to 
the IFW Image Problem Mailbox. 



